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摘 要 : 潜在 蒸 散 发 在 区 域 水 量 平衡 \ 干 旱 程 度 评价 农作物 需 水 量 等 方面 的 研究 中 具有 重要 的 作用 。 然 而 ,潜在 
蒸 散发 的 空间 化 处 理 一 直 以 来 都 是 相关 研究 面临 的 一 个 挑战 。 基 于 新 疆 地 区 1960 一 2017 年 66 个 气象 站 的 观测 数 
据 , 通 过 一 种 简单 的 参数 方程 ,实现 了 潜在 蒸 散 发 (&7, ) 的 空间 化 处 理 。 研 究 结果 表明 :人 简化 参数 方程 中 的 2 个 


重要 参数 a,c 的 空间 分 布 呈 现 一 定 的 规律 ,参数 a 在 空间 上 呈现 东南 高 ,西北 低 的 特点 ;参数 则 随 海拔 的 增加 而 


WK, © 简化 参数 方程 与 Penman-Monteith 方法 相 比 , 拟 合 结 果 在 日 .月 和 季节 尺度 上 的 R (SK 0.90, H RP 
值 随时 间 尺 度 的 增加 而 增 大 。@) 将 简化 参数 方程 的 拟 合 结果 与 CRU 数据 和 MOD16A2 数据 进行 对 比 发 现 :简化 参 
数 方程 与 Penman-Monteith 方法 拟 合 的 R 值 较 高 , 拟 合 效果 和 偏差 指标 表现 更 佳 ,而 CRU、MOD16A2 数据 拟 合 的 


R^ 值 较 低 。 简 化 参数 方程 所 获得 的 潜在 蒸 散发 精度 高 
区 潜在 蒸 散发 估算 的 简便 有 效 的 方法 。 


空间 分 辨 率 更 高 (500 m x500 m) ,是 一 种 适用 于 新 疆 地 


关键 词 : 潜在 蒸 散发 (ETo ) ; 参数 方程 ; Penman-Monteith 公式 ; 空间 插值 


洪 在 蒸 散发 是 指 具备 充分 土壤 水 分 条 件 下 热 面 
的 全 部 莹 发 量 〇 0 。 它 是 区 域 水 分 循环 和 能 量 平衡 
的 重要 环节 ,对 地 区 干旱 程度 评价 .农业 灌溉 管理 、 
水 文 模型 建立 .气候 变化 影响 等 方面 的 研究 具有 重 
BER CO), RTE AR HL Ae tL BKK VE SS TE WR 
发 ,用 ET, 表示 。ET, 的 计算 方法 较 多 ,但 认可 度 高 
且 常 用 的 方法 为 Penman-Monteith 公式 (PM ZS 
式 )59 , 它 是 联合 国 粮农 组 织 计 算 潜 在 燕 散 发 的 首 
HARO, PM 公式 基于 物理 常量 来 计算 ET, ,需要 
测量 某 地 区 的 众多 气候 要 素 , 包 括 气 温 净 辐射 日 
照 时 数 、 风 速 .空气 相对 湿度 ”] 。PM 公式 计算 的 
ET, 常常 作为 与 其 他 模型 结果 比较 的 基准 中 -2 。 但 
是 ,由 于 PM 公式 的 计算 需要 较 多 气象 观测 数据 ， 
此 ,在 一 些 缺少 资料 的 地 区 应 用 就 受到 了 极 大 的 限 
制品 。 针 对 这 种 情况 ,一 些 奉 代 性 公式 开始 出 现 。 
这 些 公式 可 概括 为 2 种 类 型 ,一 是 基于 温度 的 计算 
方法 "~; 二 是 基于 温度 和 辐射 的 计算 方法 。 
前 者 使 用 气温 数据 计算 ,后 者 则 是 使 用 气温 和 总 辐 
射 计算 ET, 值 。 然 而 ,由 于 这 些 计算 方法 发 展 于 特 
定 的 气候 条 件 和 地 区 ,其 模拟 结果 会 因 不 同 地 区 的 
气候 条 件 而 产生 偏差 ,因此 ,更 多 的 研究 开始 针对 本 
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地 ET, 模型 的 校准 工作 。Allen 等 "首次 将 这 类 模 
型 用 于 当地 ET, 估算 ,他 们 利用 当地 关于 风速 的 函 
数 校准 并 验证 了 Blaney-Criddle 模型 对 彭 曼 数 据 的 
ASME JE ALF H 2822 BU t A A HUE. 2S 
似 的 模型 还 有 Blaney-Criddle 模型 (3) | Priestley-Tay- 
lor 模型 5 Thornthwaite 模型 $ Hargreaves 参数 
Jr fU? 。 近 年 来 ,经 PM 公式 改进 的 三 参数 和 两 参 
数 模型 “…'” 因 其 计算 方式 简便 ,数据 容易 获取 , 故 
被 广泛 应 用 于 气象 站 点 数据 不 完善 或 气象 数据 获取 
难 的 地 区 的 ET, 估算 。 

除了 模型 模拟 和 校准 ET, SUE PR LE ZS AL AZ 8 
间 数 据 产 品 ,例如 GLDAS?? CRU? MOD16A2'” 
等 是 目前 使 用 较 多 的 可 靠 空间 数据 集 。 但 由 于 这 些 
数据 产品 空间 分 辨 率 和 时 间 分 辨 率 较 低 '” ,一 定 程 
度 上 限制 了 它们 的 广泛 应 用 。 针 对 此 ,有 研究 利用 
简化 参数 模型 ,通过 实现 温度 和 总 辐射 空间 化 ,进而 
对 ET, 进行 空间 插值 运算 ,得 到 了 更 高 分 辨 率 的 潜 
在 藻 散 发 的 空间 化 产品 中 。 

本 研究 旨 在 利用 简化 参数 方程 进行 新 疆 地 区 
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ET, 的 计算 与 空间 插值 ,并 通过 与 PM 公式 计算 结 
果 的 对 比 , 评 价 简化 参数 方程 在 站 点 斥 度 上 的 估算 
精度 。 利 用 决定 系数 (R ) .平均 绝对 误差 (MAE ) 、 
相对 偏差 (BIAS ) 等 指标 评估 基于 简化 参数 方程 的 
ET, 空间 数据 与 CRU 和 MODI6A2 数据 产品 的 拟 合 
精度 与 偏差 ,评估 简化 参数 方程 在 新 疆 地 区 的 适 
用 性 。 


1 数据 来 源 与 研究 方法 


数据 来 源 
本 研究 采用 新 疆 地 区 66 个 气象 站 点 (图 1) 所 
记录 的 每 日 平均 气压 .气温 、 相 对 湿度 风向 风速 和 
日 照 时 数 等 数据 。 数 据 来 源 于 中 国 气 象 数据 共享 网 
(http ://data. cma. cn) 。 新 疆 地 区 的 洪 在 蒸 散发 空 
间 数 据 分 别 来 自 时 序 为 1961 一 2014 年 的 CRU TS 
3.2 月 尺度 ET, 数据 (http ;// data. ceda. ac. uk/badc/ 
cru/ data/cru, ts/cru, ts, 3. 23/data/) 和 MOD16A2 1 
供 的 2000—2014 年 全 球 ET, 月 尺度 数据 (http :// 
files. ntsg. umt. edu/data/NTSG _ Products/MOD16/ 
MOD16A2_ MONTHLY. MERRA_GMAO_ IkmALB/ 


1.1 


干旱 区 


GEOTIFF _0. Sdegree/) , 它们 的 空间 分 辩 率 均 为 
0.5? x0.5?, 

1.2 研究 方法 

1.2.1 简化 参数 方程 及 拟 合 效果 评估 方法 ”在 气 
象 要 素 缺少 的 地 区 ,PM 公式 的 应 用 受到 限制 。 针 
对 此 ,有 研究 者 提出 仅 依 赖 于 温度 和 总 辐射 计算 
ET, 的 新 方法 9] 。 简 化 参数 方程 [以 下 简称 参数 方 
程 , 公 式 (1) ] AS PM 公式 有 一 定 关系 ,其 中 
aR, 代表 PM 公式 中 辐射 项 ,表示 输入 辐射 与 输出 
辐射 的 差 值 。(1 - cT) 是 PM 公式 中 分 母 项 的 近似 
值 "”。 参 数 方程 公式 如 下 : 
aR 
1 oP (1) 
SUP PET 指 参数 方程 计算 的 潜在 蒸 散发 值 ;了 指 该 
站 点 的 每 日 平均 气温 (% ) ;R, 是 指 理想 大 气 睛 天 下 
各 站 点 每 日 大 气 上 界 的 总 辐射 量 (MJ +m? . do'), 
与 当地 纬度 和 太阳 高 度 角 有 关 , 用 公式 (2) 计算 
得 到 。 


PET = 


sO Gd [osin(g)sin(8) +cos(@) 
cos(8) sin( c, ) | (2) 


510 680 km 


图 1 新 疆 地 区 66 个 站 点 空间 分 布 示意 图 


Spatial distribution of 66 meteorological stations in Xinjiang 
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式 中 :Gy 是 太阳 常数 等 于 82 kJ - m ?^;d, 是 日 地 平 
均 距 离 的 倒数 ;w, 为 当地 太阳 高 度 角 ;ep 是 地 理 纬 
度 ;6 为 太阳 赤 纬 角 , 以 上 3 个 变量 单位 为 弧度 
(rad) ,d, 和 6 由 日 序数 计算 得 到 ,而 co, 由 地 理 纬度 
和 日 序数 计算 得 出 LR, 的 计算 在 MATLAB 软件 
中 实现 。 参 数 akg + kJ) 和 参数 cCC 7) 是 由 联 
合 国 粮农 组 织 提 供 的 FAO — 56 Penman-Monteith 公 
式 " 计算 的 ET, 拟 合 得 到 ,采用 SigmaPlot11.0 中 非 
线性 方程 拟 合 实 现 。 公 式 如 下 : 
fs) =T e 

UP ix, 表示 每 日 总 辐射 R,;x, 表示 每 日 平均 气温 
T(?C) if(x, x, ) z& PM 公式 计算 得 到 的 每 日 ET o 
分 别 给 参数 ac 赋值 为 0. 000 1 和 0. 01 ,得 到 与 站 
点 一 一 对 应 的 参数 a scs 

PM 公式 计算 得 出 的 ET, 值 简写 为 ET, pmo 
数 方程 计算 的 ET, CET, ) 进行 对 比 ,使 用 决定 系 
BCR?) .平均 绝对 误差 (MAE) 、 相 对 偏差 (BIAS) 来 
评价 和 校准 参数 方程 在 模拟 估计 ET, 时 的 拟 合 效 
RO. Fern R? 越 接近 1 ,代表 拟 合 效果 越 佳 ,平均 
绝对 误差 MAE( mm) 与 相对 偏差 BIAS(% ) 越 小 , 表 
IR ET, 和 ET wa 间 的 偏差 值 越 小 。 
1.2.2 空间 插值 方法 ”参数 方程 实现 ET, 的 空间 
化 需要 以 下 几 个 步 又。 首先 ,使 用 Anusplin 软件 对 
每 日 平均 气温 进行 空间 插值 ,该 软件 采用 的 插值 方 
法 是 薄板 样 条 函数 法 。 然 后 ,基于 日 序数 和 地 理 
纬度 ,利用 公式 (2) 在 ArcGIS 软件 中 实现 每 日 总 辐 
HFR, 的 空间 栅 格 运算 。 

最 为 关键 的 2 个 参数 ac 空间 化 的 具体 步 又 
H:O aye 分 别 与 站 点 经 度 ,海拔 和 温度 进行 回 


(3) 


归 分 析 ,建立 参数 we 的 拟 合 方程 ;@ 基于 该 拟 合 
方程 和 相关 因子 进行 参数 a A 的 栅 格 运算 ,标记 
为 ww 和 ec;@ 以 a' 和 wc 作为 协同 变量 ,采用 样 条 曲 
线 法 进行 参数 a c 空间 插值 。 最 后 ,基于 每 日 平均 
气温 了 .总 辐射 R, 和 参数 ac 的 空间 插值 结果 , 利 
用 参数 方程 (1) HEIT ET, ,空间 栅 格 运算 ,得 到 新 
性 全 区 的 ET, 空间 分 布 。 


2 研究 结果 


2.1 不 同时 间 尺 度 下 五 Tu wo 与 ET, 的 分 析 对 比 

由 图 2 可 以 看 出 ,采用 参数 方程 计算 的 ET, moa 
与 PM 公式 计算 的 ET, ,结果 趋势 一 致 ,精度 较 高 。 
在 日 时 间 尺 度 上 ,E7, ,表现 为 高 估 的 时 间 段 主要 
为 年 初 (1 ~86 d) .年 中 (167 ~230 d) 和 年 末 (300 ~ 
365 d) ,平均 高 估 幅 度 为 0.28 mm; 其 余 时 间 段 
ET, wu 则 表现 为 低估 状态 ,平均 低估 幅度 为 0.25 
mm, JA E, ET) .表现 为 高 估 状 态 ,平均 高 估 
0.08 mm ( K| 2a) Æ H IE NR EE, ET ma 与 
ET, u 相 比 , 在 4 月 .5 月 .8 月 .9 月 10 月 表现 为 低 
fh ,平均 低估 7.47 mm ,低估 幅度 最 大 为 10.19 mm, 
与 此 相对 ,在 其 余 月 份 ETu ,ua 则 表现 为 高 佑 状态 OE 
均 高 估 7.74 mm, 最 大 高 估 幅 度 12.19 mm( 图 2b) 。 
在 季节 尺度 上 ,ET ,ww 与 £76,, 相 比 趋势 最 为 接近 ， 
总 体 上 FT， 高 估 了 5.31 mm( 图 2c)。 平均 绝对 
偏差 MAE 的 年 内 变化 (图 2a) 和 月 尺度 变化 (图 
2b) 都 随 着 时 间 的 增加 而 上 下 波动 ,1 月 4 月 7 月 、 
9 月 和 12 月 在 波峰 ,3 月 .6 月 .8 月 、10 月 处 在 波 
谷 ;并 且 季 节 尺 度 上 (图 2c),MAE 冬季 最 高 ,夏秋 
季节 最 低 。 


此 外 ,基于 站 点 数据 计算 的 多 年 日 月、 季节 的 
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图 2 不 同时 间 尺 度 PM 公式 计算 的 潜在 蒸 散发 (ET ) 与 参数 方程 7, (B CET, moa) 的 年 内 平均 变化 和 


平均 绝对 偏差 (MAE) 


Fig.2 Compared results of the potential evapotranspiration (ET, ,,,) calculated by the PM formula and the ET, value of the 
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ET, a ET, 之 间 的 决定 系数 R 评价 结果 表明 ， 
R^ 值 在 日 尺度 .月 尺度 和 季节 尺度 分 别 为 0. 991, 
0.991 .0. 993 , 拟 合 精度 会 随时 间 尺 度 增 加 而 增加 。 
因 在 月 尺度 下 ET, 4,5 ET, 在 66 个 站 点 上 的 决 
定 系 数 R 平均 值 为 0.94, 有 33. 3396 的 站 点 ET, pm 
5 ET, moat R^ 值 大 于 0.95, 有 57.09% 的 站 点 R 
值 介 于 0.92 ~0.95 之 间 , 仅 有 3.03% 的 站 点 RR 值 
小 于 0.90。 因 此 ,笔者 采用 月 尺度 下 的 7, HERE 
准 评价 参数 模型 下 的 ET, ,i 值 ,并 计算 各 个 站 点 参 
Zt a c 的 值 。 

ET, 5 ET, 在 各 个 站 点 上 的 平均 绝对 偏差 
(MAE) 与 相对 偏差 (BIAS ) 可 以 反映 出 参数 方程 与 
PM 公式 的 偏差 程度 (图 3)。 评 价 结果 表明 :) 新 
JE 66 个 站 点 的 MAE 均值 为 0.38 mm, Hif MAE 大 
于 0.50 的 站 点 占 总 数 的 12. 1296 ; MAE 介 于 0.40 ~ 
0.50 的 站 点 占 总 数 18.18% ;MAE 介 于 0.20 ~ 0. 40 
的 站 点 占 总 数 65. 15% ;而 MAE 介 于 0. 10 ~0.20 
的 站 点 占 总 数 的 4.55% 。 总 体 上 ,超过 59. 98% 的 
站 点 MAE 值 小 于 0. 40 ,因此 ,MAE 总 体 较 小 。@) 
66 个 站 点 的 BIAS 均值 为 0.016 0 ,其 中 大 于 0. 04 
的 站 点 占 总 数 1. 52% ; 介 于 0. 02 ~0. 04 之 间 的 站 
点 占 总 数 的 37.88% ;0. 01 ~ 0.02 之 间 的 站 点 占 总 
数 的 37.88% ;小 于 0.01 的 站 点 占 总 数 的 22.72% 。 
总 体 上 ,有 超过 50. 00% 的 站 点 相对 偏差 值 小 于 
0. 02 ,偏差 较 小 ,说 明 拟 合 效果 较 优 。@) zs MAE 
空间 分 布 (图 3a) 表 明 , 北 疆 大 部 分 站 点 MAE 偏 高 ， 
南 疆 站 点 MAE 相对 较 低 ,其 中 塔 中 站 MAE 最 低 ,为 
0.16 mm; 4 JE BIAS 空间 分 布 和 MAE 类 似 ( 图 
3b) , 北 疆 多 数 站 点 BIAS {r , PA Sit 2 BOW A BIAS 


平均 绝对 偏差 (MAE)/mm 
* 0.10 ~ 0.20 
e 0.20 ~ 0.30 
® 0.30 ~ 0.40 
9 0.40 ~ 0.50 
@ 0.50 ~ 0.60 


图 3 


Fig.3 Spatial distribution of mean absolute error (MAE) and relative bias ( BIAS) of 66 meteorological stations 


新 疆 66 个 站 点 在 月 尺度 下 平均 绝对 
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Ur 究 
值 偏 低 。 
2.2 参数 ac 和 ET, 的 空间 分 布 特征 

将 参数 方程 与 PM 公式 计算 的 每 月 ET, 均值 进 
行 拟 合 ,得 到 各 个 站 点 的 参数 ac 值 ,参数 a 最 大 值 
为 0. 000 1 ,最 小 值 为 0. 000 047 ;参数 e 最 大 值 为 
0. 04 ,最 小 值 为 0. 015。 运 用 样 条 曲线 法 对 2 FS 
数 进行 空间 插值 ,得 到 参数 ce 的 空间 分 布 规律 (图 
4a.4b) 。 结 果 表 明 :参数 a 的 空间 分 布 与 太阳 辐射 
相关 ,在 低 纬 度 地 区 太阳 辐射 越 高 ,a 值 则 越 大 ,而 
高 纬度 地 区 太阳 辐射 较 低 , 参 数 a 也 随 着 辐射 的 降 
低 而 减 小 ,其 空间 变化 呈现 东南 高 ,西北 低 。 而 参数 
c 则 与 海拔 相关 性 高 ,c 的 分 布 会 随 着 海拔 高 度 的 变 
化 而 变化 ,海拔 越 高 值 越 大 ,海拔 越 低 c 值 越 小 。 

由 参数 方程 空间 插值 得 到 的 新 疆 地 区 全 年 和 7 
H ET, 空间 分 布 表 明 : 

(1) 1960—2017 年 新 疆 地 区 多 年 平均 ET, 为 
1 028 mm ,并且 存 在 显著 空间 差异 (图 4c)。 从 全 疆 
来 看 ,平原 区 年 均 ET, 为 1 130 mm, 山 区 为 864 
mm ,平原 区 的 蒸发 量 远 远大 于 山区 。 平 原 区 的 全 年 
温度 高 ,降水 稀少 ,蒸发 过 程 强烈 ,而 山区 温度 随 着 
海拔 的 升 高 而 降低 ,温度 整体 较 低 , 且 降 水 较 多 3S 

(2) 南北 疆 的 多 年 平均 ET, 也 存在 明显 差异 
(图 4c) , 南 疆 地 区 年 均 ET, 值 为 1 047 mm, JU S8 Ht 
区 为 996 mm。 在 南 疆 平原 区 ,E76 达到 1 187 mm; 
南 疆 山区 潜在 燕 散 发 较 低 , 为 899 mm。 北 性 平原 区 
多 年 平均 ET, 为 1 076 mm, 山 区 为 780 mm。 新 疆 
地 区 ET, 空间 变化 趋势 差异 显著 ,主要 原因 是 : 南 吐 
地 处 暖 温带 ,有 塔克拉玛干 沙漠 ,气候 干燥 少雨 ,全 
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图 4 参数 cc 和 新 疆 地 区 年 ETu 和 7 月 ET, 的 空间 分 布 


Fig.4 Spatial distribution of the parameters a and c, and the annual ET, and ET, in July in Xinjiang 
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发 较 低 。 

(3) 蒸 散发 最 强烈 的 月 份 7 月 ,新 疆 ET, 均值 
为 164 mm, 其 空间 分 布 与 多 年 平均 ET, ais 
出 平原 区 高 ,山区 低 的 特点 E EC rp i CE 
发 量 相近 。 相 比较 全 年 ET, 的 空间 变化 ,7 月 im 
ET, 空间 变化 更 加 明显 ,分 异 也 更 显著 (图 4d) 。 
2.3 ”基于 参数 方程 的 B7\ 空间 插值 数据 的 精度 分 析 

以 PM 公式 计算 的 Er, 值 为 标准 ,分 析 了 
MOD16A2 和 CRU 数据 对 月 平均 潜在 蒸 散发 的 拟 合 
精度 (图 5$) 。 结 果 表 明 :MOD16A2 数据 的 ET, 产品 
与 ET) m AY Be BBR? 均值 为 0. 87, 由 于 
MOD16A2 数据 在 裸 地 , 芒 江 及 高 山 积 雪 层 的 缺失 ， 
仅 能 提取 46 个 站 点 的 ET, 值 ,其 中 有 31 个 站 点 RR 


s 63 个 站 点 的 R^ 值 大 于 0.90。 

T, 与 ET nog .MODI6A2 数据 和 CRU 数据 的 
R*( 表 1) 分析 结果 表明 :CRU 数据 的 c 
均值 在 三 者 中 最 高 ,为 0.95 ,其 中 有 95.45% 的 站 点 
RC 值 大 于 0. 90; 参 数 方程 仅 次 于 CRU 数据 ,R 均 
值 为 0.94, 但 有 96. 97% 的 站 R^ 值 大 于 0. 90; 而 
MODI6A2 的 R^ 均值 为 0.87, 只 有 67. 39% 站 点 R? 
值 大 于 0.90。 男 外 ,不 同 计算 方式 下 的 平均 绝对 偏差 
MAE 值 在 不 同 范围 的 站 点 个 数 也 各 不 相同 ( 表 2 ) ， 


表 1 ABET, 计算 方式 R* 值 的 分 段 站 点 个 数 
Tab.1 Number of the segmented meteorological stations 
for R^ of different ET, calculation methods 


ET, 站 点 不 同 R? 值 对 应 的 站 点 数 
计算 方式 ”总 数 <0.50 0.50 ~0.70 0.70 ~0.90 0.90 -1.00 


参数 方程 66 0 0 2 64 
值 大 于 0.90。 相 较 MOD16A2 数据 而 言 ,CRU 数据 MOD16A2 46 1 7 7 31 
产品 表现 更 佳 ,决定 系数 R 均值 为 0.95 ,66 个 站 点 CRU — 66 0 1 2 63 


决定 系数 决定 系数 
- 0.00 ~ 0.40 * 0.68 ~ 0.70 
。 0.41~ 0.81 * 0.71 ~ 0.80 
e 0.82 ~ 0.95 © 0.81 ~ 0.95 
€ 0.96 ~ 1.00 € 0.96 ~ 1.00 


(a) PM-MOD16A2 (b) PM-CRU 


图 5 PM 公式 计算 的 £7 与 MOD16A2 和 CRU 提取 的 ET, 间 决 定 系 数 的 对 比 
Fig.5 Compared results of the determination coefficients between ET, calculated by Penman-Monteith formula 


and ET, extracted by MOD16A2 and CRU 


表 2 不 同 ET, 计算 方式 MAE 值 的 分 段 站 点 个 数 
Tab.2 Number of the segmented meteorological stations 


with MAE values of different ET, calculation methods 


“4 MAE x0. 50 时 ,CRU fil MODI6A2 数据 的 站 点 数 
分 别 为 41 和 0 个 ,而 参数 方程 在 该 范围 内 的 站 点 数 
有 58 个 , 占 总 数 的 87.88% 。 各 种 计算 方式 的 相对 


Im (BIAS) 也 存在 差异 , 参数 方程 的 BIAS 为 ET, WS 不 同 MAE 值 对 应 的 站 点 数 
0.015 4, MOD16A2 数据 为 0. 21, CRU 数据 为 计算 方式 总 数 <0.50 0.50~1.00 1.00 -2.002.00 ~5. 50 

-0.006 5, 参数 方程 66 58 8 0 0 

2000—2014 年 ,CRU 数据 、MOD16A2 数据 和 参 MOD16A2 46 0 2 22 22 

CRU 66 41 17 5 3 


数 方程 与 PM 公式 计算 的 87, 月 均值 散 点 图 (图 6) 
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图 6 CRU .MOD16A2 和 参数 方程 与 PM 公式 计算 的 ET, 对 比 散 点 图 
Fig.6 Comparison scatter plot between CRU ,MOD16A2 and parameter formula with the ET, calculated by Penman-Monteith formula 


表明 :CRU 数据 和 参数 方程 拟 合 效果 明显 优 于 
MOD16A2 数据 ,而 MOD16A2 数据 出 现 了 对 ET, 值 
的 高 估 。 通 过 以 上 分 析 得 出 ,参数 方程 相 比 其 他 2 
种 数据 产品 (CRU .MOD16A2 ) ,决定 系数 R 值 大 于 
0.90 的 站 点 个 数 以 及 绝对 平均 偏差 MAE 值 小 于 
0.50 的 站 点 数 都 较 高 ,因而 拟 合 效 果 较 好 。 


3 ”讨论 与 结论 


3.1 讨论 
目前 ,潜在 燕 散发 估算 方法 主要 分 为 3 种 类 
型 :GD 基于 辐射 的 方法 ,包括 Priestley-Taylor 、 


Makkink 及 Ture 1X "77? , (9) 基于 温度 的 方法 ， 
包括 Thomthwaite 和 Blaney-Criddle ARY, (3) 
基于 辐射 ,温度 和 空气 动力 学 项 的 Penman-Monteith 
JrikU? Penman-Monteith 方法 具有 很 好 的 物理 机 
制 ” ,计算 精度 高 , 常 被 作为 ET, 估算 的 标准 。 然 
而 ,Penman-Monteith 方法 的 最 大 缺点 在 于 过 度 依 赖 
观测 数据 ,在 缺少 资料 地 区 难以 推广 和 应 用 号 。 另 
一 方面 ,区 域 ET, 的 高 精度 空间 插值 一 直 以 来 都 是 
挑战 ,目前 多 数 研 究 都 是 针对 已 知 站 点 的 插值 运 
BO A Mevicar 等 中 运用 PM 公式 估算 了 黄 
河流 域 的 ET, ,并 进行 了 空间 栅 格 运算 ,得 到 了 黄河 
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流域 月 尺度 上 ET, 的 空间 分 布 。 虽 然 一 些 数据 产品 
(如 MOD16A2 ,CRU) 可 提供 空间 化 ET, ,但 其 精度 
与 空间 分 辩 率 较 低 2 ”21 。 本 研究 所 应 用 的 参数 化 
方案 ,结构 简单 ,数据 易 获 取 ( 仅 需要 温度 和 总 辐射 
数据 ) 。 而 且 参 数 化 方程 计算 精度 较 高 ,日 月 、 季 
PREFA R 值 均 大 于 0.90。 马 亮 等 3 在 新 疆 阿 
克 苏 地 区 运用 3 种 潜在 蒸 散发 模型 计算 的 ET, 与 
PM 公式 计算 的 ET, 的 相关 系数 低 于 0.70。 相 比 
较 , 本 文中 的 研究 方法 拟 合 精度 较 高 。 除 了 拟 合 结 
果 十 分 接近 PM 公式 外 ,参数 方程 还 可 实现 ET, 的 
空间 栅 格 运算 ,与 PM 公式 的 拟 合 精度 较 高 ,而且 分 
辩 率 为 500 m ,高 于 CRU 数据 与 MOD16A2 数据 的 
空间 分 辨 率 。 利 用 参数 方程 得 到 新 疆 地 区 的 潜在 蒸 
散发 空间 分 布 ,分 布 特点 是 :平原 区 大 于 山区 ; P S 
高 于 北 疆 ,符合 新 疆 地 区 ET, 空间 分 布 特征 ,与 张 山 
清 等 1 研究 结果 一 致 。 

参数 方程 是 由 彭 曼 公式 演化 而 来 ,简化 参数 方 
程 中 参数 we 的 空间 变化 有 一 定 的 规律 。 参 数 a 和 
太阳 辐射 相关 ,呈现 出 东南 高 .西北 低 的 特点 ,与 新 
疆 总 辐射 的 空间 分 布 相 近 "” ;参数 c 与 海拔 有 关 ， 
海拔 越 高 ,参数 c 越 大 。 这 是 由 于 参数 c 与 温度 呈 
负 相 关 ,而 温度 与 海拔 呈 负 相关 ,海拔 越 高 温度 越 
低 , 因 而 参数 。 与 海拔 高 度 存 在 正 相 关 关 系 , 参 数 
a.c 在 空间 分 布 上 总 体 符合 气象 因素 的 变化 规律 ， 
具有 一 定 的 参考 价值 *。 

参数 方程 可 为 ET, 空间 插值 运算 提供 可 靠 支 
撑 , 当 某 研究 区 域 缺 乏 PM 公式 所 需 数据 集 时 ,这 种 
方程 将 有 可 能 成 为 替代 PM 公式 的 方法 "1 ,经 过 校 
准 后 的 模型 参数 更 符合 当地 的 气候 条 件 , 并 为 一 些 
表征 干旱 程度 的 指标 提供 参考 ,如 干燥 度 指数 、 
水 文 模型 等 。 综 上 所 述 ,参数 方程 的 优点 是 结构 简 
单 , 数 据 易 获 取 , 计 算 精度 高 。 在 观测 手段 不 断 进 
步 .观测 网 络 逐 渐 完善 的 背景 下 ,简化 参数 方程 在 站 
ARI DCN. BER ET, 计算 上 具有 更 大 的 应 用 价值 和 
潜力 。 
3.2 结论 

本 文通 过 新 疆 66 个 站 点 1960 一 2017 年 的 气象 
观测 数据 ,评价 和 校准 了 简化 参数 方程 的 适用 性 , 重 
点 分 析 了 参数 方程 .CRU 数据 和 MOD16A2 数据 与 
PM 公式 间 的 拟 合 程 度 , 得 出 以 下 结论 : 

(1) 参数 方程 与 PM 公式 在 日 .月 及 季节 尺度 
上 的 拟 合 效果 都 表现 较 好 ,及 值 均 大 于 0.90, 且 RR 


值 随时 间 尺 度 的 增加 而 增 大 。ZET ,与 ET, moa AI 
MAE 和 BIAS 均值 分 别 为 0.38 和 0.016 0 ,总 体 偏 
差 较 小 , 拟 合 结果 较 优 。 

(2) 参数 a.c 经 空间 插值 后 ,在 空间 分 布 上 呈 
现 一 定 规律 ,参数 a 呈现 东南 高 ,西北 低 的 特点 ,与 
太阳 辐射 的 空间 变化 规律 相似 , 随 着 太阳 辐射 的 降 
低 而 减 小 ;参数 e 的 空间 分 布 则 随 着 海拔 高 度 的 增 
加 而 增 大 。 通 过 参数 方程 的 空间 插值 得 到 新 疆 ET 
的 空间 分 布 总 体 表现 为 南 疆 大 于 北 疆 , 平 原 区 大 于 
山区 的 格局 。 

(3) 评价 了 参数 方程 .CRU 数据 和 MOD16A2 
数据 与 PM 公式 的 拟 合 效 果 及 偏差 程度 ,参数 方程 
在 拟 合 效 果 和 空间 分 辨 率 方面 表现 更 佳 ,R 大 于 
0. 90 ,空间 分 辨 率 为 500 m x 500 m, 可 以 提高 ET, 
的 计算 精度 和 空间 分 辩 率 。 参 数 方程 适用 于 新 疆 地 
区 的 潜在 蒸 散 发 计算 ,为 气象 资料 缺乏 地 区 提供 了 
一 种 结构 简单 ,数据 容易 获取 的 ET, 估算 新 方法 。 
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Abstract; Potential evapotranspiration plays an important role in estimating the regional water balance and crop 
water demand and also assessing drought. However, the spatial interpolation of potential evapotranspiration is always 
a challenge for relevant researches. In this study ,a simple parametric model was used to achieve the spatialization 
of the potential evapotranspiration (ET, ) in Xinjiang based on the observed data from 66 meteorological stations in 
Xinjiang from 1960 to 2017. The results indicated that :GD There was a certain regularity of spatial distribution of the 
two important parameters a and c of the simplified parametric model. The parameter a was spatially high in the 
southeast but low in the northwest of Xinjiang. The parameter c, however, was increased with the increase of alti- 
tude ;(2) Compared with the Penman-Monteith method ,the R? values of the simplified parametric model at the daily , 
monthly and seasonal scales were all higher than 0. 90 , and they were increased with the increase of time scale; (3) 
The compared results between the fitting results of the simplified parametric model and the CRU data and 
MOD16A2 data revealed that the R? values of the simplified parametric model and of the Penman-Monteith method 
were high ,the fitting results and errors were better, but the R^ values fitted with CRU and MODI6A2 were low. Mo- 
reover ,the accuracy of potential evapotranspiration and the spatial resolution (500 m x 500 m) obtained by the sim- 
plified parametric model were high. Therefore , the simplified parametric model is a convenient and effective model 
for calculating the ET, at regional scale , and can be applied to estimate the potential evapotranspiration in Xinjiang. 
Key words: potential evapotranspiration (ET, ) ; parametric equation; Penman-Monteith method; spatial interpo- 
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